Transport properties of graphene across strain-induced nonuniform velocity profiles 
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We consider the effect of uniaxial strain on ballistic transport in graphene, across single and 
multiple tunneling barriers. Specifically, we show that applied strain not only shifts the position of 
the Dirac points in reciprocal space, but also induces a deformation of the Dirac cones, and that 
both effects are of the same order on the applied strain intensity. We therefore study the deviations 
thereby induced on the angular dependence of the tunneling transmission across a single barrier, as 
well as on the conductivity and Fano factor across a single barrier and a superstructure of several, 
periodically repeated, such sharp barriers. Our model is generalized to the case of nonuniform 
barriers, where either the strain or the gate potential profiles may depend continuously on position. 
This should afford a more accurate description of realistic 'origami' nanodevices based on graphene, 
where 'foldings' are expected to involve several lattice spacings. 

PACS numbers: 73.20.Mf, 62.20.-x, 81.05.ue 



I. INTRODUCTION 

Graphene is an atomically thin, two-dimensional layer 
of carbon atoms arranged according to a honeycomb 
lattice. After having being speculated since long as 
the ideal building block of graphite and other sp^ car- 
bon compounds, it has been recently obtained in the 
laboratory^, thereby kindling an extraordinary outburst 
of experimental as well as theoretical research activity^''^. 
Reduced dimensionality and its peculiar structure con- 
spire towards the formation of low-energy quasiparti- 
cles, which can be described as massless Dirac fermions 
with a cone dispersion relation in reciprocal space around 
the so-called Dirac points K, K', and a linearly vanish- 
ing density of states (DOS) at the Fermi level. This 
is reflected in several electronic properties already in 
the non-interacting limit, e.g. Klein tunneling^^~, the 
reflectivity^, the optical conductivityiSrii^ and the plas- 
mon dispersion relationi^"— . 

Graphene is also remarkable for its exceptional me- 
chanical properties, as is generic for most carbon com- 
pounds. For instance, notwithstanding its reduced di- 
mensionality, graphene is characterized by a relatively 
large tensile strength and stiffness^-, with graphene 
sheets being capable to sustain elastic deformations as 
large as 20%22r^. Larger strains would then induce 
a semimetal-to-semiconductor transition, with the open- 
ing of an energy gap2^"— , and it has been demonstrated 
that such an effect critically depends on the direction of 
applied straini^i^. The effect of uniaxial strain on the 
linear response electronic properties of graphene has been 
studied on quite general grounds^. 

Recently, it has been suggested that graphene-based 
electronic devices might be designed by suitably tailor- 
ing the electronic structure of a graphene sheet under 



applied strain'^^. Indeed, a considerable amount of work 
has been devoted to the study of the transport proper- 
ties in graphene across strain-induced single and multiple 
barrier&2^^. There, the main effect of strain has usually 
been considered to be that of shifting the position of the 
Dirac points in reciprocal space. However, it has been 
demonstrated that a nonuniform space variation of the 
underlying gate potential would result in a modulation 
of the Fermi velocity^i^i^. 

Here, we show that both effects are of the same or- 
der on the applied strain intensity, and should there- 
fore be considered on the same ground, when studying 
the transport properties of strained graphene. We shall 
therefore explicitly consider not only the strain-induced 
displacement of the Dirac points in reciprocal space, but 
also a strain-induced deformation of the Dirac cones, re- 
sulting in a strain-dependent anisotropic Fermi velocity. 
Specifically, we will consider tunneling through a single 
strain-induced sharp barrier, possibly subjected to a gate 
potential, and through a superstructure made of several 
such barriers, periodically repeated. More interestingly, 
we will generalize our results to the problem of trans- 
port through a tunneling structure, characterized by a 
nonuniform variation of both the Fermi velocity and of 
the gate potential, as can e.g. be brought about by a 
continuous deformation or applied uniaxial strain. 

The paper is organized as follows. After introducing 
our model in Sec. [Hi we discuss the effect of a strain- 
induced modulation of the Fermi velocity on the angular 
dependence of the transmission across a single sharp bar- 
rier, as well as on the conductivity and Fano factor for 
ballistic transport (Sec. We then consider the case 
of several such barriers, arranged in a periodic fashion 
(Sec. lIVI) . In Sec.|Vl we generalize our results to the case 
of nonuniform strain across a smooth barrier. Finally, 
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in Sec. I VII we summarize and give directions for future 
investigation. 

II. MODEL 

In unstrained graphene, low-energy quasiparticles can 
be described by the linear Hamiltonian in momentum 
space 

= hv^la ■ p, (1) 

where vp is the Fermi velocity, cr ~ (cti, (T2), with ai and 
Ti {i = 1,2,3) Pauli matrices and I the identity matrix 
associated with the two-dimensional spaces of the sub- 
lattices {A and B, say), and of the two valleys around 
the Dirac points (K and K'), respectively. Eq. ([TJ acts 
on the four-component spinors^i^i 

*p - (*A,K(p),^'i3,K(p),^'B,K'(p),-*A,K'(p))^, (2) 

where p is measured from the Dirac point one is referring 
to. Here and below, a superscript zero denotes absence of 
strain. The effect of uniaxial strain in real space is that 
of modifying the lattice vectors as 8i = {1 + e) ■ {£ = 
1,2,3), where 6^°^ = a{VS,l)/2, = o(-y3, l)/2, 

5^^^ = o(0, —1) are the relaxed (unstrained) vectors con- 
necting two nearest-neighbor (NN) carbon sites, with 
a = 1.42 A, the equilibrium C-C distance in a graphene 
sheelfS., and e is the strain tensor^- 

e=^e[{l-,y)l + {l + iy)Am, (3) 

where 

A{e) = cos(26')cr3 + sin(26l)CT^, (4) 



where the Pauli matrices now are understood to act on 
vectors of the two-dimensional direct or reciprocal lattice. 
In Eq. (|3]), 9 denotes the angle along which the strain is 
applied, with respect to the x axis in the lattice coordi- 
nate system, e is the strain modulus, and v is Poisson's 
ratio. While in the hydrostatic limit — ~1 and e = el, 
in the case of graphene one has v = 0.14, as determined 
from ab initio calculations^S,^ to be compared with the 
known experimental value u — 0.165 for graphite^S. The 
special values 9 = and 9 = tt/2 refer to strain along the 
zig zag and armchair directions, respectively. 



The possibility of describing the effects of strain 
through Eq. ([3]), i.e. elastically, implies that applied 
strain does not induce any irreversible process or me- 
chanical failure of the graphene sheet, such as disloca- 
tions, grain boundaries, or cracks. In fact, such dra- 
matic effects are not expected for strain below 20 %, 
as is predicted by calculations within density functional 
theory^^'**', and confirmed experimentally by means of 
atomic force microscopy (AFM)''^. 



In momentum space, the effect of uniaxial strain on 
the Hamiltonian Eq. ([T]) is likewise accounted for by the 
strain tensor, Eq. This is usually described as a shift 
in momentum space of the location of the Dirac points. 
However, starting from the more general, tight-binding 
Hamiltonian^, expanding to first order in the strain mod- 
ulus, and to second order in the impulses, one may show 
that applied strain also induces a deformation of the 
Dirac cones, at the same (first) order in e. Explicitly, 
one finds 
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(5) 



where acq = {a./2t)\dt/ da\ k, 1.6 is related to the loga- 
rithmic derivative of the nearest-neighbor hopping t at 
e = 0. 

Our model is based on the tight-binding approximation 
for the band structure, including only nearest-neighbor 
hopping. To this level of approximation, one does not ob- 
serve any strain-induced modification of the work func- 
tion. In order to include also such effects, one needs 



to consider also next-nearest neighbor hopping^. Mak- 
ing use of the expression for the hopping function be- 
tween two neighboring carbon p-orbitals involved in a 
TT bond, as a function of the bond length £, VppTr{t) = 
^gg-3.37(^/a-i)^ with to = -2.7 eV^, one finds 



?=%/3a 



1.7 eV X e, (6) 
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viz. a scalar term, going linear with the strain modulus e, 
whose order of magnitude agrees with the ab initio results 
of Ref. [l^. At any rate, the work function, Eq. (jG]), can 
be absorbed in an effective scalar potential U, which we 
conventionally refer to as to a gate potential below. 

Another effect that is not explicitly considered in our 
model is the deformation of the tt orbitals due to off- 
plane bending, as would be e.g. generated by an AFM 
tip. However, a change in the hopping parameters due to 
the bending of the graphene sheet can be described as an 
effective in-plane strain'*^. Specifically, one may expect 
that strain induced by an AFM tip would be character- 
ized by cilindrical symmetry, which is beyond the scope 
of the present work, where only linear barriers are con- 
sidered. We note in passing that other efficient ways to 
realize controllable strain consists in depositing graphene 
on top of deformable substrates^i^. 

The spectrum of the strained Hamiltonian, Eq. ([S]), 
is still linear, but now around the shifted Dirac points 
qdO = ±(Koe(l + ly) cos(26'), -ko£(1 + v) sin(26'))^. To 
first order in the wavevector displacement q = p T Qd 
from such shifted Dirac points, one finds 



H = hv-pcr 



q', 



(7) 



where 



[{I- Ke{l~v))l~ Ke{l + v)A{e)]ci, (8) 



and K 



Ko — ^. However, it is convenient to work in 
the reference frame with the x axis along the direction of 
applied strain. This is accomplished by a rotation in the 
sublattice AB space, described by the unitary matrix 



1 

e-'' 



(9) 



so that 



H = nvFU\9) [(7i(l - A,e)fe + f72(l - Xye)qy] U{e), 

(10) 

where A^: = 2k, Aj, = —2kv. After the rotation, Eq. 
the location of the Dirac points is given by 

qdo = ±(Ko£(l + i')cos(36i),-Koe(l + ;^)sin(36'))^. (11) 

The density operator can be expressed as 

p(r) = *^(r)*(r), (12) 

where ^'(r) = (27r)~^ J dP'ke^'^^'^'^ili). Correspondingly, 
the current density operator can be derived as^ J = 
-f [-f^^r], yielding 

J,(r) = -ez;F*t(r)(l - \e)U\e)a,U{e)-^{v). (13) 

In the following, for the sake of definitiveness, we shall 
restrict to the valley K only, thus having q = p — qo • 
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FIG. 1: One-dimensional single tunneling barrier along the x 
direction. Region II (0 < a; < D) is characterized by applied 
strain e along the 9 direction, as well as by a gate voltage Vg. 



III. TUNNELING ACROSS A SINGLE 
BARRIER 



Potential barriers for single quasiparticle tunneling in 
graphene are conventionally designed by suitably chang- 
ing the underlying gate voltage. Recently, it has been 
suggested that an equivalent effect may be induced by 
local uniaxial strain^id^. Therefore, we start by consid- 
ering a strain-induced one-dimensional step-like barrier, 
characterized by uniaxial strain applied along the direc- 
tion 9, with respect to the x axis, Eq. ([3]), with strain 
modulus s for < a; < D, and zero otherwise. Corre- 
spondingly, the Hamiltonian and current density vector 
are given by Eqs. (fTO|) and (fT3|) . respectively. In addition, 
for the sake of generality, we may also consider a nonzero 
gate potential Vg within the barrier (Fig.[T]). 

Since we are interested in stationary solutions and the 
strain-barrier is uniform along the y direction, the energy 
E and the component ky of the wavevector of an incoming 
wave are conserved. We look therefore for solutions of the 
stationary Dirac equation of the form 



( W {e)ilji{x)e'''yy , x<0, 
tpix, y) = I U'' {e)i^ii{x)e^^yy , < x < D, 
[uH9)TPm{x)e'''yy , x>D, 



(14) 



where 
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In Eqs. dH]) 
respect to the barrier 

{\E\lhv^)%\i\ip, [E 



(/9 denotes the angle of incidence with 

fcr = (|i?|/fi,UF) cos fc„ — 



fi2„2[(l _ A,e)2g2 + (1 
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9 = e = 7t/2 e = o e = 7t/2 
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FIG. 2: (Color online) Dependence on the incidence angle 
of the tunneling transmission T, Eq. p7|) . Left panel refers 
to strain applied along the zig zag direction {9 = 0), and 
(a) £ = 0.03, Ug = meV; (b) e = 0.03, Ug = -20 meV (the 
strain-induced deformation of the Dirac cone is neglected); (c) 
£ — 0.03, Ug = —20 meV. Right panel refers to strain applied 
along the armchair direction (6 — 7r/2), and (a) e — 0.01, 
Ug = meV; (b) e = 0.01, Ug — meV (the strain-induced 
deformation of the Dirac cone is neglected); (c) e — 0.01, 
Ug = -20 meV. 

\yef{ky - goy)^], s' = sgn{E - Ug), with Ug = -eVg 
Propagating waves correspond to real values of qx, while 
evanescent waves correspond to having purely immag- 
inary. 

Given the stationary character of the solution, the con- 
tinuity equation implies that V • J = everywhere. In 
particular, (J) = (?/;|J|?/;) may only depend on x, there- 
fore (Jx) is constant. The latter condition implies, at the 
barrier boundaries, 

MO') = (1 - A,£)-i/Vii(0+Xl6a) 
{1-Xxe)-'^^M{D-) = Viii(I?+). (16b) 

Enforcing the above conditions in Eqs. ([T5|). one eventu- 
ally finds for the tunneling transmission, T ~ \t\'^, 

J. ^ cos^ f 

cos^ ip cos'^{qxD) + (1 — ss' S sin (p^ s'm^ (qxD) ' 

(17) 

where qy = ky~qr,y, qx = {I- Xx£) ^\{E - Ug)"^ /fi^vl - 
(1 - A^e)2g2|i/2^ C (1 - Xxe)hvFqx/\E - Ug\, S = 
(1 - Xye)hvFqy/\E -Ug\. 

A. Angular dependence 

In order to discuss the dependence of the tunneling 
transmission on the incidence angle ip, we preliminarly 
observe that propagation within the barrier is allowed 
whenever 

?l\|(l - XyS)^{ky - q^iyf <{E- Ugf, (l8) 
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FIG. 3: (Color online) Same as Fig. g] but for E = 150 meV 
and D = 100 nm. 

where ky = (E/hv-p) simp. Within such a range, one has 
moreover total transmission (T = 1) whenever 

qxD = rnr, (19) 

n being an integer. Eq. (fTS)) differs from the usual condi- 
tion for propagation across strain-induced barriers'^^ in 
that we are not only considering a shift of the Dirac 
point qd, but also a strain- induced deformation of the 
Dirac cone, here exemplified by the substitution vf i— > 
i;f(1 - Xys). 

Figs. [2] and [3] show our results for the tunneling trans- 
mission T = T{tp), Eq. p7|) . as a function of the incidence 
angle tp, for £; = 80 meV, D = 100 nm (Fig. ^ and 
E = 150 meV, D = 100 nm (Fig.E]). In both figures, left 
(resp., right) panel refers to uniaxial strain applied along 
the zig zag {9 = 0; resp., armchair, 6 = 7r/2) direction. 

In the case of strain applied along the zig zag direc- 
tion {9 = 0, Figs. [2]and|3l left panels), curves (b) neglect 
a strain-induced deformation of the Dirac cone. Com- 
parison with curves (c), where such a deformation is 
fully included, shows that the effect of a strain-induced 
anisotropy of the Fermi velocity is that of shifting the 
angular location of the maxima {T ~ 1, Eq. p^ ) of 
the tunneling transmission. Such an effect becomes more 
important with increasing energy (from Fig. [5] to Fig. [3]), 
while the number of peaks increases, Eq. (|19l) . and the 
angular range in which the propagating regime is allowed 
widens. The effect of a strain-induced deformation of the 
Dirac cone is even more dramatic in the absence of a gate 
potential [Ug = meV, curve (a)]. Indeed, in such a 
case, neglecting the Fermi velocity anisotropy for strain 
applied along the zig zag direction would yield a uniform 
tunneling transmission T = 1, for all incidence angles ip, 
whereas we find that transmission via propagating waves 
is allowed only for \(p\ < arcsin[(l — Xye)~^], with small 
oscillations below T = 1 within, and evanescent waves 
beyond that range. A similar analysis applies to the case 
of strain applied along the armchair direction (9 = 7r/2, 
Figs. [5] and |31 right panels), which is characterized by 
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FIG. 4: Schematic top view of a graphene layer contacted by 
metallic leads, as considered in Sec. IIII Bl 



an asymmetric transmission T = T{ip), with pronounced 
oscillations for > close to the propagating edge. 

The origin of such an asymmetry of the iy9- dependence 
of the transmission can be traced back to the particu- 
lar Dirac cone vertex qd, whose shift is here considered. 
Global symmetry would be restored upon inclusion of 
the other Dirac cone. In that case, one would obtain the 
same picture, but with ip i— > ^ip. It should be emphasized 
that the stationariety condition, Eq. (fTO]). characterizes 
the occurrence of peaks in the transmission T(ip) in any 
case. In addition, for a potential barrier, in the absence 
of strain, one also recovers complete transmission (T = 1) 
at if = (Klein tunneling). 

Summarizing, at variance with previous studies^, 
from Eq. ([T7)) one obtains that the overall effect of a 
strain-induced deformation of the Dirac cones is that of 
shifting the transmission peaks, and of reducing the range 
in if at which transmission takes place. 



B. Ballistic transport 

We now consider a more realistic device, viz. a 
graphene strip of length L and width W, subjected to 
two leads at a distance D (Fig. |4])^i2iii&ii2. Following 
Ref. [3, we assume that W/L ^ 1, and that the gate 
potential within the strip is much less than the potential 
of the leads, |V^| <C |Vl|- Moreover, we assume that the 
graphene strip be characterized by uniaxial strain, with 
modulus £ and strain direction 9, and explicitly consider 
the deformation of the Dirac cones induced by the ap- 
plied strain. The energy levels of Ref. \^ are therefore 



modified into 

E = Ul 



shvY^kl + kl, a;<0, x>D, (20a) 



<x < D, 



(20b) 



where again % = ky — q^y, and Ul = —e-Vl and Ug = 
—eVg. The limit \Vl\ — > oo is equivalent to the limit 
<y9 — > 0, and the transmission, Eq. (|17p . reduces to 



1 



cos'^iqxD) + f]{ky) sin^ (q^D) ' 
for propagating waves in the valley a = K, where 

(E-Ug)^ 



Ijiky) = 



[E - Ug)^ - h^V^il - Xye)^ky i (/D^ ) ^ 



(21) 



(22) 



and the minus {resp., plus) sign applies to the valley 
a = K {resp., a = K'). Analogous expressions hold 
for the transmission T^^'^{ky) in the evanescent case, 
with ri{ky) > —ri{ky), cos{qxD) i— >■ cosh.{qxD), and 
sm{q,j.D) 1-^ sm]i{qxD). The transmission for a general 
(propagating or evanescent) wave therefore reads 

T^ky) = e[viky)]Tr''iky) + (1 - eivikyWr^'iky), 

(23) 

where Q(t) is the Heaviside (step) function. Integrating 
over ky and summing over both valleys, one obtains the 
conductance across the barrier (Landauer formula) "^^i*^^ 



G = 



2e 



27r 



Ta{ky), 



(24) 



where the factor of 2 takes into account for the spin de- 
generacy, the conductivity 



and the Fano factor—!^ 



F = 1 



Z^Q J-oo 271- ^a(l^y) 



(25) 



(26) 



In Eq. (j25p for the conductivity, the summation over 
the valleys contributes with an additional factor of two, 
whereas this factor cancels in the definition of the Fano 
factor, Eq. CT . 

Before discussing our results, let us observe that the 
inclusion of a strain-induced deformation of the Dirac 
cone in the expressions of the conductivity, Eq. (j25p . and 
of the Fano factor, Eq. ([25]) . amounts to the replacements 



D ^ Deff 

E ^ E,s -: 



CE, 



(27a) 
(27b) 
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for the strip width and incident energy, respectively, in 
the corresponding expressions, cr^°^ and F^^"^, say, with- 
out cone deformation, with 



1 - A^e' 
1 

1 — XyE 



(28a) 
(28b) 

(29) 

As a consequence, while ImiE-i-ou'^'^^D^E) — ie^/nh, a 
universal constant*^, in the presence of applied uniaxial 
strain one finds 



In particular, one explicitly finds 



1 

lim a{D,E) = -^. 



(30) 



Only in the case of hydrostatic strain {ly — ~1, — Xy, 
^ = 1) does one recover the universal limit, regardless 
of the strain modulus^. On the other hand, one finds 
liniE ->-o F{D, E) = i, corresponding to strongly sub- 
Poissonian noise^, regardless of applied strain. 

In the opposite limit, the conductivity across a sin- 
gle barrier in the absence of strain is linear in energy, 
cr(°) Si {e^ /h)D\E\/hvF for E ^ oo, with damped oscil- 
lations characterized by a pseudoperiod AE such that^ 
DAE/hvp = TT. In the presence of strain, such results 
arc modified by Eqs. ([SS)) . so that a{E) « (Toa{E) for 
E oo, with 



4e2 D\E\ 



-c, 



(31) 



with damped oscillations characterized by a pseudope- 
riod given by 



Tivy 



(32) 



In view of the fact that \Xx\ > |Aj,|, one may conclude 
that applied strain induces a slight change in the slope 
of (T vs while it modifies the pseudoperiod of the 
oscillations more substantially. 

Fig. [5] shows our results for the scaled conductivity 
in the presence of uniaxial strain (e — 0.03 — 0.15) ap- 
plied along the armchair direction [6 = 7r/2). When 
the conductivity (j{E) is normalized with respect to its 
asymptotic limit, Eq. ([3T|) . and plotted against energy 
E scaled with the strain-dependent pseudoperiod AE, 
Eq. ([32|) . results corresponding to different values of the 
strain modulus collapse into a single curve, displaying 
damped oscillations, as prescribed by Eq. ((32|) . Simi- 
larly, Fig. |6] reports our results for the Fano factor as a 
function of scaled energy. Again, the results for all the 
strain moduh here considered (e — 0.03 — 0.15) collapse 
into a single, oscillating curve. Note that the universal 
limits F{E = 0) = ^ and F^o = Htxie-^oo F{E) = \ are 
recovered in all cases, regardless of applied strain. Such 
results do not depend on the direction 9 of applied strain. 
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FIG. 5: Conductivity across a graphene strip {D = 100 nm) 
normalized to asymptotic large-energy behavior, Eq. (|3ip . vs. 
energy scaled to the pseudoperiod, Eq. p2|) . Actually shown 
are four curves, all collapsing into a single one, corresponding 
to strain applied along the armchair direction (6 = 7r/2), with 
£ = 0.03, 0.05, 0.10, 0.15. 




FIG. 6: Fano factor for ballistic transport across a graphene 
strip. All parameters are as in Fig. [5] Dashed lines represent 
the universal low- and large-energy asymptotic values, F{0) — 
4 and -Foo = |, respectively. 



IV. TRANSMISSION ACROSS MULTIPLE 
BARRIERS 



We next consider quasiparticle tunneling across N 
identical barriers, each of width £, two nearest neighbor 
(NN) barriers being separated by the distance £, such 
that 2N£ = D (Fig. [7]). We assume a position-dependent 
strain modulus e{x) and gate potential energy U{x), with 



e{x) = e_, {m — \)l < X < ml, (33a) 
= £+, m£ < x < (m -f 1)£, (33b) 
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where r, t {resp., r' , t') are the amplitudes of the reflected 
and transmitted waves in region I (resp., Ill), cf. Fig. [T] 
Indeed, one explicitly finds^i^ 



(it)- 

-{tr 



r'it') 

{t'Y 



(40) 



Therefore, for the conductance across a single barrier, 
one finds 



D 



FIG. 7: Schematic plot of the multiple barrier, as considered 
in Sec.HVl 



and 



U{x) = (to - 1)£ < X < m£, (34a) 

= U+, mi<x<{m+ l)e, (34b) 

with TO = 1, . . . 2N — 1. We further consider the possibil- 
ity of contacting the two extrema of the chain of barriers 
with leads at the potential Vl- Eqs. ([T5)) then suggest to 
look for a solution of the Dirac equation in the form 



i;{x,y)^U\0) 



0(x) 



a/1 - Xxe{x) 



Jkyy 



(35) 



so that (/'(x) is a continuous function at the barriers' 
edges. The stationary Dirac equation for 4'{x) can then 
be casted in the form of an evolution equation^, so that 
(t>{x) = T^-^"^ {x,xq)(/){xo), where the evolution matrix 
T*^^^ (a;, a;o) in turn obeys the equation 



j^TrW(^,^o) = 

ax 



iq^^le{x)T^l + 



i E-U{x) 



hvp 1 — Xxs{x) 



1 - Xxe{x) 



tW(x,xo), (36) 



withTW(a;o,xo) =1. For a single barrier, the evolution 
matrix is related to the transfer matrix by^ 



l^'\x,xo) - g;i(^)TW(x,xo)Q.(¥'), (37) 



where 



V2 Ve 



1 



(38) 



includes the incidence angle tp of the incoming spinor, 
Eq. (|15ap . and s — sgn (E). In the limit of metallic leads 
(|Vl| co), one has 0, with Q+(0) = ^(f^^ +<7x), 

Q+\0) = Q+(0), and Q_(0) = Q+(0)a,, Qz\0) = 
axQ+{0)- The elements of the transfer matrix can be fur- 
thermore related to the elements of the scattering matrix 
across the barrier, 



S = 



r t' 
t r' 



(39) 



G 



2e2 
h 



Tr (th) = ^Tr ((M^^^^^W^) , (41) 



where Tr = W^^l°^^dkyl2-K. Correspondingly, the 
transmission for an incoming quasiparticle with energy 
E and transverse wavevector ky in valley a is Ta{ky) = 

(M^i^%iW\ and the expressions for the conductivity, 
Eq. (pS]) . and Fano factor, Eq. (^5)) . follow straightfor- 
wardly. 

The solution of Eq. (l36l) for the transfer matrix is de- 
rived analytically in Appendix both for a single and 
for a multiple barrier, in presence of strain-induced de- 
formation of the Dirac cone. Making use of Eqs. (IA10[) 
for the transmission Ta{ky) in Landauer's formula for the 
conductivity, Eq. (|25|) . and in the definition for the Fano 
factor, Eq. (pS)) . one again finds that the conductivity 
in strained graphene, and strained graphene where the 
strain-induced velocity anisotropy has been neglected, 
are related by means of Eqs. ([28|) . ([29| . but now with 
D = 2Ni, and 



c 



2(e+ + e-), 

^(C+ + C-), 

1 - \ye± 
1 - Xxe± ' 
1 

1 - Xye± ' 



(42a) 
(42b) 
(42c) 
(42d) 



Eq. (|30l) in the hmit E 
wardly, with ^ given now by Eq 



then follows straightfor- 



Moreover, the 
conductivity at large energies is characterized by an over- 
all linear behavior, interrupted by dips with decreasing 
depth, which result from the coherent superposition of 
the damped oscillations produced by scattering off the 
edges of the single barriers. The energies En at which 
such dips occur are asymptotically given by (cf. Ap- 
pendix |A| 



(43) 



with n an integer. 

Fig. [5] shows our numerical results for the conductivity 
in strained graphene, with strain applied nonuniformly 
along the armchair direction, across a superlattice of 

= 10 barriers. At variance with Fig. [51 we have not 
scaled a with its asymptotic behavior at large energies. 
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E/E, 




FIG. 8: (Color online) Conductivity a-{E) in units of uq — 
vs. energy E, scaled with respect to the approximate 



location of the first dip, Ei, as given by Eq. (|43|) . Subsequent 
dips then occur close to integer values of the ratio E/Ei. 
Uniaxial strain is applied along the armchair direction {9 — 
it/2) in the case of a multibarrier superlattice, with = 10 
barriers, £ = 25 ran {D = 500 nm). Different curves refer to 
nonuniform strain moduli within and outside NN barriers (cf. 
Fig.0, with (a) e+ = 0.004, e_ = 0; (b) e+ = 0.003, £_ = 0; 
(c) e+ = 0.002, e- = -0.001; (d) e+ = 0.002, e_ = 0.001; (e) 
e+ = 0.0005, e_ = 0. In all cases, we set U± — 0, for the sake 
of simplicity. 



FIG. 10: (Color online) Conductivity o-{E) in units of ao = 
Ae^/h, vs. energy E, scaled with respect to Ei, as given by 
Eq. (|43[l . Uniaxial strain is applied along the zig zag direction 
(9 — 0) in the case of a multibarrier superlattice, with A'^ = 10 
barriers, £ = 25 nm {D = 500 nm). Different curves refer to 
nonuniform strain moduli within and outside NN barriers (cf. 
Fig.0, with (a) e+ = 0, £_ = 0; (b) e+ = 0.03, £_ = 0; (c) 
£+ = 0.05, e_ = 0; (d) e+ = 0.07, £_ = 0; (e) e+ = 0.10, 
£_ — 0. In all cases, we set U± — 0, for the sake of simplicity. 
Inset shows the conductivity scaled with respect to its large- 
energy asymptotic limit, a/aoo, as a function of scaled energy, 
E/Ei. 




I ' ' ' ' ' ' ' ' 1 

01 23456789 

E/Ei 

FIG. 9: (Color online) Fano factor F vs. scaled energy 
E/Ei, for transport across a multibarrier superlattice, with 
nonuniform uniaxial strain applied along the armchair direc- 
tion {9 = 7r/2). All parameters are as in Fig. [S] Inset shows 
the universal low-energy asymptotic behavior in the various 
cases. In the limit E ^ 0, the universal asymptotic value, 
F{0) — ^, is recovered. 



Eq. (|31|) . As expected, the overall linear behavior of (j{E) 
is interrupted by dips, whose approximate energy loca- 
tion is given by Eq. (|43| . While such dips get damped as 
energy increases, they are nonetheless enhanced with re- 
spect to the case in which the strain-induced deformation 
of the Dirac cones is neglected^, especially those corre- 



sponding to even integer values of n in Eq. (I43p . Corre- 
spondingly, the Fano factor (Fig. [5]) is characterized by 
essentially analogous features, with bumps occurring at 
approximately £"„, Eq. (|43p . In particular, the univer- 
sal limit at low energy, F{0) — i, is recovered as in the 
single-barrier case, regardless of applied strain. 

Fig. [To] shows our numerical results for the conductiv- 
ity in strained graphene, but now for nonuniform strain 
applied along the zig zag direction. At variance with the 
armchair case (Fig.[8l), for strain applied along the zig zag 
direction the conductivity seems not to be characterized 
by prominent dips as a function of energy. This may ex- 
plained by a reduced coherent superposition of the effects 
due to each single barrier. However, if the trailing lin- 
ear dependence on energy is divided out (Fig.[Tni inset), 
one may again recognize 'oscillations', with extrema ap- 
proximatively occurring at En, as given by Eq. (|43p . At 
variance with the armchair case, the Fano factor exhibits 
a strain-dependent asymptotic limit, for large energies 
(Fig. Ilip . with increasing deviations from the unstrained 
behavior Foo = |, with increasing strain modulus e (at 
least within the strain range that has been numerically 
investigated). On the other hand, both the oscillations 
as a function of scaled energy E/Ei and the low-energy 
limit F{Q) = i (Fig. [TTJ inset) are recovered. 
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FIG. 11: (Color online) Fano factor F vs. scaled energy 
E/Ei, for transport across a multibarrier superlattice, with 
nonuniform uniaxial strain applied along the zig zag direc- 
tion {9 = 0). All parameters are as in Fig. 1101 Note the 
deviations from the large-energy asymptotic limit for the un- 
strained case, = I (dashed line). The low-energy univer- 



sal limit, F(0) 
less of strain. 



(inset, dashed line), is recovered, regard- 




FIG. 12: Schematic single tunneling barrier, with smooth 
strain profile, Eq. (|48|) . Dashed line depicts a sharp barrier, 
corresponding to the limit a 0. 



V. TRASMISSION ACROSS A SMOOTH 
BARRIER: EFFECT OF CONTINUOUS STRAIN 



Although considerable insight is afforded by analyti- 
cal solutions to the problem of tunneling across single or 
multiple sharp barriers, there is sufficient evidence, both 
experimental^^ and theoretical"^^, that barrier edge ef- 
fects are also important to determine the transport prop- 
erties across corrugated graphene. Here, we therefore 
consider the case in which uniaxial strain is applied in a 
nonuniform but continuous fashion to a graphene sheet, 
which can be modeled by a single barrier with smooth 
strain and gate potential profiles, e = e{x) and U = U{x), 



respectively. Such a description includes and generalizes, 
in particular, a continuous Fermi wavevector profile, as 
considered in Ref. [s^ . 

On quite general grounds, one may expect that a 
smooth potential profile (whether induced by strain or 
by gating) introduces a new length scale, a say [as in 
Eq. P5)l below], which is the linear size over which the 
potential strain varies appreciably. Such a new length 
scale has then to be compared with the atomic scale, 
measured by the lattice step o, on one hand, and with 
the Fermi wavelength Ap = hvp / {27rE) corresponding to 
the incident energy E, on the other. The approximation 
of a sharp barrier (no smoothing) then holds whenever 
a <C a <C Ap, i.e. at sufficiently large incident ener- 
gies. On the other hand, the detailed structure of the 
barrier needs to be considered when a ~ Ap. In both 
cases, we are interested to the more general and realis- 
tic cases where a a, where one may additionally ne- 
glect the occurrence of K-K' coupling. Indeed, truly 
sharp electrostatic barriers on the order of the electron 
wavelength are quite difficult to be realized, as is e.g. 
demonstrated by the occurrence of Fabry-Perot oscilla- 
tions of the conductance in graphene heterostructures as 
narrow as ~ 20 nm, where a resonant cavity is formed 
between two electrostatically created bipolar junctions^. 
Such oscillations are more accurately described when the 
smooth structure of these potential barriers is taken into 
account, whereas intervalley scattering can be safely ne- 
glected (see Supplementary Information in Ref. 55). An- 
other instance of nonuniform barrier, where smoothing 
effects are important, is the strain-induced ripples super- 
lattice experimentally realized in Ref. \^ which smooth- 
ing is essential on a length scale of ~ 100 nm, whereas 
intervalley processes are negligible. 

The kinetic part of the Hamiltonian for graphene sub- 
jected to uniform strain e along the direction 9 is 



H = u\e)a,hv, ( iv, 



(44) 



where Vi = wf(1 — A^e), and summation over the re- 
peated index i — 1,2 is understood. In order to gen- 
eralize Eq. (|44l) to the case of a nonuniform, but con- 
tinuous strain profile e = e{x), one may be tempted to 
perform the replacements Vi i-> Vi{r) = vp[l — Xie{x)] 
and Qd qD(r), Eq. ([TT|). with e ~ e{x). However, the 
resulting Hamiltonian must be symmetrized, in order to 
preserve hermiticity, thus leading to the model Hamilto- 
nian for a nonuniform strain profile: 



H = U^0)a, 



1 



hvi{r) -Vi - gDi(r) 



U{0). 



Vi - qoiir) hvi{v) 



(45) 



Eq. (j45p includes the effect of nonuniform, continuous 
strain both as a shift in the position of the Dirac points, 
and as a deformation of the Dirac cones (nonuniform 
and anisotropic Fermi velocity), at variance e.g. with 
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Ref. I35I where a nonuniform velocity is considered, but 
an isotropic profile is assumed. As in the case of a single, 
sharp barrier fSec. lIII| . continuity of the current density, 
Eqs. (IT51) . suggests to seek for a solution of the stationary 
Dirac equation in a form analogous to Eq. ((35|) . viz. 



JkyV 



One explicitly finds [cf. Eq. (IMl) ] 

d4>{x) 



(46) 



dx 



1 - Xyejx) / _ (0) . ^^ 
1 - X^e(x) T'' l^y'^"^)) 



+i 



e{x) 
E-U{x) 
(1 - \x£{x)) hvF 



(x). (47) 



We have solved Eq. (I47p numerically, for the nonuni- 
form, smooth strain profile 



e{x) 



£0 



1 



1 



tanh(D/4a) VI + e^^^/'^ 1 + e-(^--D)/'» 



as shown in Fig.[T3 Such a strain profile is essentially fiat 
for \x—D/2\ <^ a, where e{x) « Eq, and for \x—D/2\ :s> a, 
where s{x) w 0. In the hmit a/D 0, Eq. (gS]) tends 
to the sharp barrier considered in Sec. IIIII Therefore, 
asymptotically for |a;| — )■ 00, the solutions of Eq. (jTZj) 
must merge into Eqs. (|15p . in regions I and III. We have 
therefore taken an initial value (f)(x — xq) in the form 
of Eq. ()15c|) . for xq — 5D, and integrated Eq. (jTZ)) back- 
wards for a; ^ 0. Comparing the numerical solution with 
Eq. (|15a|) . one may extract the refiection coefficient r, rel- 
ative to an incident wave with unit amplitude incoming 
from a; > 0, as the Fourier weight with respect to its 
negative frequency component, whence the transmission 
T{(j)) follows straightforwardly. As a cross-check of our 
procedure, we have also verified that the continuity equa- 
tion, Eq. (|16p. holds true, within the numerical error. 

Figs. [13] and [14] show our numerical results for the tun- 
neling transmission T{ip) across the smooth strain bar- 
rier, Eq. (|48| . with = 100 nm and different values 
of the smoothing parameter, a/D. Fig. [13] refers to an 
incidence energy E = 80 meV, corresponding to an inci- 
dent wavelength Ap = hv-p / {2'kE) « 1.3 nm. One finds 
that transmission of propagating waves is allowed for in- 
cidence angles Lp such that i^cr- ^ ^ci-+, with 



± arcsin 



1 



1 — Xy£0 

in the zig zag case (0 = 0), and ip > ipa-, with 
1 hvp 



1 — A,,£o 



+ T^^ok(1 - I') 



(49) 



(50) 



in the armchair case {9 — 7r/2), independent of the 
smoothing parameter a/D. Outside that window, trans- 
mission takes place via evanescent waves only, and 
T{ip) sa 0. For strain applied along the zig zag direction 



e = o 



e = 7t/2 




0.5 



-0.5 




0.5 



FIG. 13: (Color online) Tunneling transmission vs inci- 
dence angle (p across a smooth strain barrier, Eq. (|48|l . with 
D — 100 nm, and incidence energy E = 80 meV (Ap — 
hv-p/{2nE) ~ 1.3 nm). Left panel refers to strain applied 
along the zig zag direction {9 — 0), with eo — 0.1. Right 
panel refers to strain applied along the armchair direction 
(48) (6 = 0), with £0 = 0.01. In both cases, the different lines 
correspond to different values of the smoothing parameter, 
VIZ. (a) a = (sharp barrier); (b) a — lO^^D = 1 nm; (c) 
a = 10~^-D — 10 nm. In all cases, U{x) = 0, for the sake of 
simplicity. 



e = o 



e = 7i/2 




FIG. 14: (Color online) Same as Fig. [T3 but with E 
150 meV (Ap f» 0.7 nm). 



(9 — 0, Fig. [131 left panel), Eq. (|49)) predicts the existence 
of critical angles I^Scril < 7r/2. This is a direct conse- 
quence of the strain-induced deformation of the Dirac 
cones [Aj, 7^ in Eq. (j^U)) ]. Both in case of strain ap- 
plied along the zig zag and armchair directions, increas- 
ing the smoothness parameter a/D away from the limit 
of a sharp barrier {a/D — 0) suppresses the oscillations 
in T{(p) within the propagating window, until a > Ap, 
in which case transmission is almost undisturbed by the 
presence of the barrier. These results are confirmed by 
Fig. 1141 where we consider quasiparticles with larger in- 
cident energy i? = 150 meV, corresponding to a smaller 
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Fermi wavelength Ap ~ 0.7 nm. While the transmission 
window widens and the number of oscillations increases, 
smoothening the strain profile immediately washes out 
the deviations of the tunneling transmission from unity. 
In ending this section, we note that the procedure applied 
to extracting the tunneling transmission from the numer- 
ical solution of Eq. (j47p can be generalized, in principle, 
to the case of an arbitrary nonuniform strain potential, 
such as a superlattice of several smooth barriers, such as 
Eq. (gSD. 
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Appendix A: Transfer matrix across a multiple 
barrier 



VI. CONCLUSIONS 



We have studied the effect of a strain-induced modu- 
lation of the Fermi velocity on several transport proper- 
ties of graphene, such as the angular dependence of the 
tunneling transmission, the conductivity, and the Fano 
factor. After considering the cases of a single sharp tun- 
neling barrier, and of a superstructure of several, peri- 
odically repeated, such sharp barriers, we have specifi- 
cally studied the case in which both the modulus of ap- 
plied uniaxial strain, and possibly an applied gate poten- 
tial, depend continuously on position. This is expected 
to afford a more accurate description of real 'origami' 
device^, in which 'foldings' of a graphene sheet would 
conceivably involve several lattice spacings. In the case 
of sharp tunneling barriers, we have demonstrated that 
the effect of a strain-induced deformation of the Dirac 
cone is of the same order of the strain-induced shift of 
the Dirac points, and should therefore be taken into ac- 
count on the same basis. In particular, we have found 
that strain modifies the quasi-period in energy that reg- 
ulates the occurrence of dips in the conductivity across 
a superstructure of several sharp barriers, due to coher- 
ent scattering off their edges. Such effect is however less 
dramatic in the energy dependence of the Fano factor. 
Finally, we have generalized our results to embrace the 
case of a generic nonuniform strain, and possibly a gate 
potential, profile. Besides allowing a more accurate anal- 
ysis of tunneling transmission across smooth barriers, es- 
pecially at low incident energies, which are expected to be 
more sensitive to local deviations from uniformity, such 
an approach can be applied to describe arbitrary strain 
superstructures, albeit numerically. 

Among the already available experimental results 
which could be described in terms of a strain-induced 
deformation of the Dirac cones, we mention Raman 
spectroscopy^ and the strain-dependence of both the 
longitudinal and the recently predicted transverse plas- 
mon mod o ^°i'~. Moreover, transmittance measurements 
with polarized light between the near-infrared and the 
ultraviolet on uniaxially strained graphene may provide 
information on the Dirac cone deformation^!^. 



In the case of a single barrier (N = 1, 2£ = D), Fig.[71), 
Eq. p6l) for the transfer matrix admits the analytical 
solution 



(AO) 



exp 



i (E ~ Ug)D 

exp I : — a,. 



+ ^^(^.-9^)^-.), (Al) 



corresponding to the initial condition M^^^ (0, 0) = I, and 
to a uniform strain e and to a gate potential energy 
Ug across the barrier. The second matrix exponential 
in Eq. (|A1|) can be made more explicit, by making use 
of the following identity for a linear combination of the 
Pauli matrices. 



exp(a • cr) 



sinh a 



a 



I cosh a, 



(A2) 



where a = {J2^ af)^^"^, and a, € C (i = 1, 2, 3). 

We next consider a single barrier, but now with nonuni- 
form strain modulus and gate potential energy, i.e. 
e{x) = e_ and U{x) = U- within the barrier {0 < x < £), 
and e{x) — £+ and U{x) — U+ beyond the barrier's sec- 
ond edge {£ < X < 2£; cf. Fig. [71). In this case, one finds 
M(i)(2^,0) = M+{£)M_{£), where M±(^) are given by 
Eq. (|Aip , with D !->■£,£ i-i' e±, and U U±. One finds 



(2£,0) 



■ (0)/ 1 



(A3) 



where Mi is a unimodular matrix, detMi = 1. Specifi- 
cally, one finds 



A + ij], 



(A4) 



where 



sinh(g_^) sinh((7+€) , , 
A = — w_it+) 

q- q+ 

+ cosh((7_£) cosh(g4.^), 
r] = — sinh((7_^) cosh{q^£) 

q- 

sinh(g_|_^) cosh(g_^), 

9-1- 



(A5a) 
(A5b) 
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with 



Since A = ^TrMi = coshz, one finds explicitly 



K± 
U± 

q± 



1 — Xye± 



1 - Ke± 
E-U± 



(A6a) 
(A6b) 
(A6c) 



whence Eq. (|^T|) follows straightforwardly. 

Finally, in the case of barriers {D = 2Ni, Fig. \7}, 
iterating Eq. (IA3|) N times, one has 



(AO) 



■ (0) , 



N 
1 I 



(A7) 



where for the A^th power of the unimodular matrix Mi 
one may use an identity due to Chebyshev^^, and specif- 
ically obtain 



11 smhz smhz 



Here, we have denoted the eigenvalues of Mi by e^^, with 
z € C Finally, one finds for the transmission 



T{ky) = [(mW(D,0))*^ (mW(D,0))^^ 



cosh^(A^z) + 



A2 - 1 



sinh^(Afz) 



.(A9) 



Ta{ky) — T^^°'^{ky) 



cos\Ny) + ^^sin^ {Ny) 



(AlOa) 



with y — arccosA, if |A| < 1, 



rr^'iky) 



cosh^(Aj;) + sinh^(7Va;) 

A^ — 1 



(AlOb) 



with X = log |A + VA2 - 1|, if |A| > 1, 



[1 + ry^A^- 



2 at2]-1 



(AlOc) 



if |A| = 1. In particular, one finds A ^ cos(u_|-^ + u^£), 
ioi E ^ oo, whence Eq. p3|) follows. 
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